Theory of dynamic crack branching in brittle materials 



E. Katzavt*, M. Adda-Bedia^ and R. Arias* 
Labor atoire de Physique Statistique de I'Ecole Normale Superieure, CNRS UMR 8550, 
24 rue Lhomond, 75231 Paris Cedex 05, France. 
^Departamento de Fisica, Universidad de Chile, CasiUa 487-3, Santiago, Chile. 

February 6, 2008 



Abstract 

The problem of dynamic symmetric branching of a tensile crack propagating in a brittle 
material is studied within Linear Elastic Fracture Mechanics theory. The Grifhth energy 
criterion and the principle of local symmetry provide necessary conditions for the onset of 
dynamic branching instability and for the subsequent paths of the branches. The theory pre- 
dicts a critical velocity for branching and a well defined shape parameterized by a branching 
angle and a curvature of the side branches. The model rests on a scenario of crack branching 
based on reasonable assumptions and on exact dynamic results for the anti-plane branching 
problem. Our results reproduce within a simplified 2D continuum mechanics approach the 
main experimental features of the branching instability of fast cracks in brittle materials. 



1 Introduction 



The continuum theory of fracture mechanics is concerned with the quantitative description of the 
mechani s ms of crack nucleatio n, the conditions under which they propagate and their dynamics 



(jFreundl . 119901 : iBrobergj . Il999l ). For brittle materials, the relationship between internal stress 
and deformation and the balance laws of physics dealing with mechanical quantities do not 
include the possibility of material separation. Consequently, the equation of motion of the crack 
front is based on additional statements on crack growth. The most frequently used criterion of 
crack propagation in two dimensional elastic brittle materials consists of two parts: Griffith's 
hypothesis and the principle of local symmetry. 
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The Griffith's energy criterion ([Griffithl . Il92d : iFreund llQQOl : iBrobera . Il999l ) states that the 
intensity of the loading necessary to produce propagation is given by ^ = F, where G is the 
energy release rate and F is the fracture energy of the material (i.e. the energy needed to create 
new surfaces) . The principle of local symmetry states that the crack advances in such a way that 
in-plane shear stresses always vanish in the viciiiity of the crack tip. This rule was first proposed 
for quasi-static cracks (iGol'dstein and Salganikl . [iQTJ; lLeblondl.ll989l). and generalized to rapidly 



moving cracks (jAdda-Bedia et al 



1999^). It was shown in (lAdda-Bedia et al 



19991 ^ that the 



two criteria arise from the same physical origin. The energy release rate is the component of the 
driving force along the direction of crack motion, Fi. Griffith's energy criterion may then be 
reinterpreted as a material force balance between Fi and a force that resists the crack advance, 
i.e. Fi = F. However, this equation is not sufficient to determine the trajectory of a crack. If one 
assumes that material force balance holds at the crack tip, one should impose the component of 
the material force perpendicular to the direction of crack propagation to vanish. This condition 
is identically satisfied if the loading in the vicinity of the crack tip is purely tensile. 



The Griffith criterion and the princi ple of local symmetry predict adequately the path and the 
stabi l ity of slowly propagating c r acks (lAdda-Bedia and Pomeaul . Il995l : lAdda-Bedia and Ben Amail . 



1996 



Bouchbinder et al 



20031 : 



confirm the theoretical results (Yuse and Sano 



Mardei 



2004') . Controlled experiinents o n quasi-static cracks 



19931 : 



Ronsin et al 



crack propagat i on, experiments on different brittle niateria. l s (IRavi-Chandar and Knauss 



Fineberg et al 



199 



199S 



1992 



Gross et al 



Livne et al 



1993 



Boudet et al. 



1996 



19951). In the case of fast 



Sharon et al 



1991 



1984 



Sharon and Finebergj . 



20051 ) have identified a dynamic instability related to a transition from 



a single crack to a branched crack configuration. The instability occurs when the crack speed 
exceeds a critical velocity Vc, which does not depend on the applied traction and on the ge- 
ometry of the plate. Above Vc, a single crack is no longer stable. Instead, a repetitive process 
of micro-branching occurs, which changes the crack dynamics : the acoustic eraission from 



the crack increases (jBoudet et al.l . 1 19961 : 



Boudet and Ciliberto 



2000; 



Gross et al 



19931 ). the 



crack velocity develops strong oscillations and a pa ttern, which is correlated with the veloc- 



ity oscillations, is observed on the fracture surface (jFineberg et al 



Sharon and Finebergl . Il996l . 



199S 



Livne et al. 



20051 '!. 



1992 



Sharon et al. 



1995 



Some aspects of this dyn a.mic instabil i ty wer e described in the framework o f the theory of 



brittle fracture mechanics ( Adda-Bedia 



2004b 



2005 



Bouchbinder et al. 



20051 ). These studies 



were based on Eshelby's approach which states that, as in the single crack case, a growth 
criterion for a branched crack must be based on the equality between the energy flux into the two 
propagating tips and the energy req uired to open the material and create new surfaces as a result 



of this propagation (jEshelbvl . Il970l ) . The problem of determining the in-plane dynam ic stress 



200 



3). 



It was 



intensity factors immediately after branching was formulated in (jAdda-Bedial . 
shown that the in-plane elastic fields immediately after branching exhibit self-similar properties, 
and that the corresponding stress intensity factors do not explicitly depend on the velocity 
of the single crack tip before branching. These properties are similar to the mode III crack 
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branch ing problem, which was solved exactly in (jAdda-Bedia and Ariasl . 120031 : lAdda-Bedial . 
2004al ) ■ This similarity suggests that under plane loading configurations, the jump in the energy 
release rate due to branching is maximized when the branches start to propagate very slowly. 
Under this assumption, the branching of a single propagating crack under tensil e loading was 



found to be energetically possible when its speed exceeds a certain critical value (jAdda-Bedia 



20051 ). The theoretical results for the c ritical velocity and the branching angle agree fair l y wel l 



Livne et al. 



2OO5I I. 



with the available experimental results (jFineberg et al.l . Il992l : ISharon and Finebergl . 



1996 



199S 



The main purpose of the present study is to perform a quantitative analysis of the subsequent 



path s followed by the branches. Following ([Karihaloo et al 



1981 



Leblondl . Il989l : lAmestov and Leblondl . 



1992) , the asymptotic expansion of the stress field at the tip of a curved extension of a branched 



crack is presented. Using these exact results, the paths selected by the br anched cracks are de- 



1996, 


1999: 


Livne et al.. 


2005) 



20051 ) is recovered without introducing any additional parameters. The 
present study shows that both the branching instability threshold, the branching angle and the 
subsequent paths of the branches can be predicted within the continuum 2D theory of brittle 
fracture mechanics. Note that the present analysis provides a necessary condition for branching 
and not an in stability mechanism for it. In addition, the branching instability in real systems is 



of 3D nature (Livne et al 



20061 ). Therefore, the present study should be seen as a step towards 



a complete understanding of this phenomenon. 



For the sake of clearness, the second section of this paper summarizes the results obtained for the 
branching problem: the branches shape, their dynamics as well as the dynamic instability are 
determined. In this paper a systematic analysis of the branching pr oblem is made, and due to 



reasons of completeness, some of the results already commun icated in (lAdda-Bedia 



2004b 



2005!) 



are presented here again. Higher order terms than those in (jAdda-Bedial . l2005l ) were calculated, 
as well as numerical corrections of lower order terms are given. The third section of the paper 
presents the detailed study of the branching problem. The fourth and last section solves an 



elastostatic problem related to the experimental setups (IFineberg et al.l . Il992l : 



1995 



Sharon and Finebera . 119961 . Il999l : 



Livne et al 



Sharon et al 



2OO5I ). where the interest is to determine 



possible outcomes of the non singular T stress at the original crack tip. 



2 Summary of results and future prospects 



First, a static analysis of crack bra nching under plane loading conditions is done. It is an 
exact approach that follows that of (lAmestoy and Leblondl . Il992l ) for the kinking case. Using 
the principle of local symmetry as a criterion for determining crack's trajectories, as well as 
asymptotic analysis of the local stress configurations at the crack tips of branches of a given 
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longitude and characteristic shape, the crack's branching angle and subsequent curved paths are 
determined. This is done by solving integral equations for the elastic potentials, with the help 
of conformal mapping and perturbation techniques. 

Secondly, the dynamic crack branching is addressed by arguing that the plane loading case should 
not differ very much qualitatively from the exact solutions of anti-plane branching, where indeed 
dynamic branching results can be safely retrieved from the static ones. Much of this is based 
on the argument that it is plausible that the branching occurs with the new branches starting 
at vanishing speeds since the elastic energy release rate at each crack tip is maximal in that 
case. Under these assumptions, application of Griffith's energy criterion leads to the critical 
speed for branching (once all the post-branching elastic quantities are replaced by their static 
counterparts). 

These results, detailed in the following sections, support the following scenario for the process 
of dynamic branching instability in brittle materials 



• The critical velocity at which the crack tip can branch depends on the material parameters 
only through the fracture energy and the elastic constants. This prediction for the critical 
velocity agrees with the available experimental results. 

• The paths that the branches take consists of two universal features: a branching angle 
of 27° followed by a curved extension described by a single curvature parameter which 
is calculated below. The branching angle as well as the general shape of the extensions 
coincide with those seen in the experiments. 

• The velocity of the branches is vanishingly small right after branching, which seems to 
be a peculiar characteristic. However, since the predicted acceleration of the tips after 
branching is very high, this may explain the experimental results. 



The above scenario for a single branching event can now be integrated into a general picture 
of the propagation dynamics of the crack, which can reproduce the fractography of the broken 
sur face, by considering in ad dition the stability analysis of the branched configuration reported 



in (Bouchbinder et al. 



20051 ). It turns out that the symmetric form of branching is unstable, in 



the sense that sooner or later one extension continues to grow while the other one slows down 
until it stops. This leads to a pattern of a broken surface composed of a sequence of branching 
events, where each time only one extension survives. The surviving crack tip accelerates between 
branching events and then decelerates abruptly in the n ext branching even t. Although the 3D 
nature of the instability is not taken into account here (|Livne et al.l . 120061 ) , we think that this 
scenario gives a coherent physical interpretation of the fractography of broken samples. In order 
to describe completely the observed patterns, our 2D analysis should be coupled to an instability 
mechanism of the crack front itself. However, we expect that features such as the critical velocity 
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for branching or the branches' shape should be modified by the 3-D nature of the problem as a 
secondary effect. 



In the following section we detail the results that led to the above mentioned conclusions. 



2.1 Stress field ahead of curved extensions of a branched crack 

Consider an elastic body containing a straight crack with symmetrically branched curved exten- 
sions of length d. and a branching angle Avr (see Fig. [1]). Let XOY denote the coordinate system 
with the OX axis directed along the initial straight crack, and let Y1OY2 denote the coordinate 
system with the OYi axis directed along the tangent to the upper extension at the point O. 
These two coordinate systems are obviously related by 



Y2 



X cos Avr + y sin Avr , 
Y cos Xt: — X sin Att . 



(1) 
(2) 



Following the approach developed in (jKarihaloo et al 



1981 



Leblond 



19891 : 



Adda-Bedia 



2004b|) 




Figure 1: Schematic representation of a straight crack with two symmetrically branched curved 
extensions. 



it can be shown that the asymptotic shape of the crack extension is necessarily given by 



Y2 = aY^I'^ + 0{Y^ 



(3) 



where a is a curvature parameter whose dimension is (length) Moreover, the expansion of 
the st atic stres s inten sity factors K[{^s) (Z = 1, 2) at the crack tip in powers of I obeys the general 



form (Leblond 



19891 ) 



= E Fim{X)Km+ E [Gm{X)T6i^ + aHUX)K„,]Vi + 0{e) . (4) 



m=l,2 



m=l,2 



In this expansion, Ki and T are the static stress intensity factors and the nonsingular stress 
in the universal expansion of the stress field at the original crack tip O without the branched 
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extensions. Ki and T are given by 



0"i 



E 

1=1,2 



'2'nr 



T.'i {e)+T6ix6,x + 0{V7) 



where (r, 6) are polar coordinates referred to the branching point O, and S 



(0 



(5) 



are known functions 



describing the angular variations of the stress field components (jBrobera . Il999l ) . The functions 
Fim^ Gi and are universal in the sense that they do not depend neither on the geometry of 
the body nor on the applied loading. They depend only on the b ranching angle and the i r corn - 
putation can be performed following the approach developed in (jAmestov and Leblondl . Il992l ). 
Note that the asymptotic expansions given by Eqs. ©-([Ij) are applicable to crack extensions 
obtained by actual propa gation of the initial crack and not simply by arbitrary machining of 
the body ( Leblondl . Il989l ). Due to the linearity of the problem, the expressions ©-(jl]) can be 



predicted from dimensional arguments. Since the Ki^s scale as stress x vlength and T scales 
as stress, the first order expansion of the stress intensity factors in dH) must involve an addi- 
tional parameter whose dimension is 1/ ^length. This parameter is provided by the asymptotic 
expansion ([3]) of the branched extension. 
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Figure 2: The elements of the matrix Fi^^X) for the symmetrically branched crack (solid lines) 
and for the kinked crack (dashed lines). 

In Figures [21 [3] and H] we present the results for the universal functions Fim, Gi and Hi^ for 
the branched cracks. As a comparison, we also present the results obtained for these functions 
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Figure 3: The elements of the vector G/(A) for the symmetrically branched crack (solid lines) 
and for the kinked crack (dashed lines). 



for the kinked crack problem (lAmestoy and Leblondl . Il992l ). Th e functions F i m for a symmet- 
rically branched configuration have been already computed in (|Adda-Bedial . l2005l ). while the 
computation of Gi and Hi^ is new. Details of this analysis are given in the following section. 
Note that th e results for the func tions Fi^ in Fig. [2] correct numerical inaccuracies in the results 
reported in (jAdda-Bedial . l2004bl ) . These inaccuracies have of course a quantitative implication 
by correcting certain values such a s the critical branch ing angle. However, more importantly, it 
contradicts the statement made in (jAdda-Bedial . l2004bl ) concerning a possible difference between 
the predictions of the principle of local symmetry and the maximal energy release rate criterion, 
as will be explained below. 



Once a detailed expansion of the stress intensity factors is available, it remains to be com- 
bined wit h a propagatio i i criteri o n, in order to get cr ack path prediction. Griffith's energy 



criterion (iGriffith 



192d: 



(jGol'dstein and Salganik 



Freund 



1990 



Brobera . Il999l ) and the principle of local symmetry 



19741 : iLeblondl . Il989l ) imply 



r 
. 



(6) 
(7) 



where fx is the Lame shear coefficient of the material. Note that Eq. ([7]) imposes the symmetry 
of the stress field in the vicinity of the crack tip which in turn affects the crack direction of 
propagation. Therefore, the crack path is mainly selected by the principle of local symmetry, 
while Eq. ^ controls the intensity of the loading necessary to advance the crack. In the 
following, the stability of a tensile crack and the path selection of branched cracks will be 
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Figure 4: The elements of the matrix HimiX) for the symmetrically branched crack (solid lines) 
and for the kinked crack (dashed lines). 



discussed in view of these general results. 



2.2 The shape of the branched cracks 



Consider a straight crack initially und er pure tensi l e (mod e I) loading that branches into 
two symmetrical cracks. As shown in (jAdda-Bedial . l2004bl ) . the paths of the branches can 
be determined if each crack ti p satisfies the principle of local symmetry during propagation 
(IGol'dstein and Salganikl . Il974l ). Imposing this (Eq. ([7|)), the following two conditions that 
determine the branching angle A and the curvature parameter a should be satisfied 



F2i{X) - 
G2{X) 



= 0, 
T 



H2i{X)Ki 

Using the results presented in Figures [21 [3] and [U Eq. ([8]) gives A = Xr 



(8) 
(9) 

0.15 corresponding to a 



Adda-Bedia 



bran ching angle of 27°, in agreement with the results of (llsida and Noguchi Il992l : 
20051 ). Also we find that G2(Ac)/i^2i(Ac) = —1.16, and consequently the sign of a is determined 



by the sign of the ratio T/Ki. Therefore, the convexity of the branches' paths depends on the 
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sign of the T-stress: If T < the branches will tend to come back towards the initial direction 
of the crack prior to branching (see Fig. [5]) , while if T > the branches will diverge away from 
this direction. 




Figure 5: Path followed by the branch when a < 0, and a comparison with the functional form 
Y{X) ~ X^/^ proposed in (jSharon and Finebergj . Il999l ). 



Before continuing, it is useful to define the material parameter k = (q/cs)^, where q and Cs are 
the dilatational and shear wave speeds. This quantity w ill serve to corn pare theoretical results 
with experiments, k is related to the Poisson ratio u by (iBrobergl . Il999l ) 



— =K 



for plane stress 
1Z2u plane strain. 



(10) 



In order to estimate the T-stress for the rea l experimental setups of (jFineberg et al.l . Il992 



Sharon and Finebergl . Il999l : 



Livne et al. 



20051 ). we solved the elastostatic problem of a semi- 



infinite straight crack in an infinite strip of half- width W, whose top and bottom edges are fixed 
at positions ±5 above their initial positions (see Fig. [T3|) . The analysis has been performed for 
both clamped edges and shear free edges boundary conditions. The details of the calculations are 
given in section It turns out that for typical situations the ratio T/Ki is negative and so the 
branches will be stable, i.e. they will not diverge away from the initial direction of the primary 
crack (see Fig. [HI and Fig. [T5b). For the typical value of k = 3, we get T/Ki = — 0.43/\/W 
when the edges of the strip are clamped, which leads to a value of the curvature parameter 
a = — 0.5/\/W- This theoret ical value is consistent with t he est imation of the parameter a from 



the experimental results of (jSharon and Finebergj . I1996I . Il999l ). However, a more quantitative 
comparison with the experiments would require to take into account finite size effects, the three 
dimensional geometry of the sample and dynamical effects. 



It is interesting to mention that when the edges are shear free, the sign of the T-stress can be 
made positive by applying an additional loading Too by stretching the strip in the X-direction. 
This results in the existence of a critical Tcr{n) such that the path of the branches becomes 
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unstable, i.e. diverges away from the direction of the primary crack (see Fig. 115b). Such a 
situation where different boundary conditions can change the stabihty of the propagating crack 
is of interest, and can certainly be tested experimentally. 



2.3 The dynamic branching instability 



Until now, only static aspects of the branching instability have been discussed. In order to 
address dynamical aspects of this problem, such as determining the onset of dynamic branching, 
let us consider the following scenario: A semi infinite straight crack that propagates at a speed 
v{t) for t < suddenly stops at t ^ 0^. At t ^ 0^, the crack branches locally with a branching 
angle equal to Att. For t > 0, the n ew branches p r opagate at a velocity v'{t) in the new 



19751; 



Freund 



1990 ) that the dynamic stress 



directions ztAvr. It is well established (jKostrovl . 
intensity factors, Ki{t,v), of the straight crack prior to branching are related to the rest stress 
intensity factors, Ki(t, 0), of the same configuration by 



Kiit,v) = kiiv)Ki{t,0), 



fill 



where ki{v) {I = 1,2,3) are known uni versal functions of the instant aneous crack tip speed v{t) 



whose explicit forms can be found in (jFreund 



( Freund . 


1990: 


Broberg. 


1999) 



Since in this problem there is no time scale, and consequently no length scale, against which 
the independent variables can be scaled, the dynamic stress intensity factors immediately after 
branching, K'i{0~^ ,v' ,v) can always be written in the form of a universal function of the veloc- 
ities and branching angle, multiplied by the static stress intensity factors immediately before 
branching, Ki{0~,v = 0), i.e. 



,v ,v) 



J2kiv')J'imi\v',v)KiiO-,0). 



fl2) 



As in the quasi-static case (|Leblondl . ll989l ^. the matrix f is universal in the sense that it depends 



neither on loading configuration nor on the geometry of the body. Indeed in the limit that is 
considered here, the dynamic branching problem does not involve radiation effects, so it is 
always equivalent to a crack propagating in an unbounded body. Moreover, T should approach 
the elastostatic solution for a vanishingly small velocity of the side-branches, namely 



lim Tim{X,v',v) 
ti'— >o 



(13) 



Based on the solution to t he anti-plane branch ing problem (jAdda-Bedia and Ariasl . l2003l : lAdda-Bedial . 
2004al ) ■ it was shown in (jAdda-Bedial . l2005l ) that the dependence of the velocity of the single 
crack tip before branching is suppressed from the stress distribution that has to be balanced 
during the propagation of the branches. Consequently, the matrix elements J^im related to plane 
loading situations should also be independent of the velocity prior to branching, namely 



(14) 
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In order to proceed, one must come up with a growth criterion for a branched crack. It is 
well estabhshed that the dyri amic energy release rate G for a single straight crack is given by 



( Kostrov 



19751 : 



Freund 



1990) 



where 



3 3 

G = -^J2Ai{v)Kf{t,v) = ^Y.9iiv)Kf{t,0) , 
1=1 1=1 



gi{v) = Aiiv)kfiv) 



(15) 



(16) 



The functions Ai{v) and gi{v) do not depend on the details of the applied loading, nor on the 
configuration of the body being analyzed. They only depend on th e local instantaneous speed 
of the crack tip and on the properties of the material (jFreundl . 119901 ). Fig. [B] shows the function 
gi (v) that will be used in the following. 



A growth criterion for a branched crack must also be based on the equality between the elastic 
energy flux into each pr opagating tip a nd the energy that is used in creating new broken surface 
during this propagation ([Griffithl . 1 1920 ) . The dynamic energy release rate is a quantity associated 
to a single moving crack tip, and so it has to be determined for each crack tip. When the primary 
crack before branching is under pure mode I loading and due to the symmetry of the branching 
configuration, the energy release rate immediately after branching G' for each crack tip is given 
by 



G' 



1 

2^ 



giiv')J^f,{X,v')+g2iv')J^i,{\v') KfiO-,0). 



(17) 




Figure 6: The universal function gi{v) for k 
wave speed. 



3. Here and elsewhere, cr denotes the Rayleigh 



According to the generalized Griffith's criterion (iGriffithl . Il920l ). the crack must grow in such 
a way that the energy release rate is always equal to the dynamic fracture energy of the ma- 
terial, r(f), which is assumed to be a property of the material and whose value may depend 



on the instantaneous crack tip speed (jFreundl . 119901 : 



Boudet et al 



1996 



Sharon and Fineberg 



11 



19991 ). This growth criterion should hold for the crack tips before and after branching, and so it 
introduces a relation between the energy release rates immediately before and after branching, 
namely 



G' 



T{v) 



G, 



(18) 



which is a nec essary condition for the existence of a branching configuration. If this condition 
is not fulfilled then single crack tip propagation would be maintained. Let us stress again that 
this condition does not provide an instability mechanism to the branching process. However, 
whichever instability mechanism, it should respect Eq. (jlSp . 



Under in-plane configuration, the exact dependence of the matrix on the crack velocity af- 
ter branching and on t he branching angle is not available. However, the exact resolution of 
the mode HI problem (jAdda-Bedial . l2004al ) does give an indication about its general behav- 
ior, since in many cases, physical aspects of crac k propagation of correspondin g anti-plane and 



in-plane config urations are qualita tively similar ([Broberg 



the results of ( Adda-Bedia 



199S 



Freund 



1990 ). In particular. 



2004al ) show that T33{X,v') depends only weakly on v': the ratio 



^ssi'^j'v')/ F3:i{X) is very close to unity (up to 5%) for all values of A and v' . We are then led 
to assume that this property will also hold for all the matrix elements J-'im{X,v'). Therefore, 
the energy release rate immediately after branching for in-plane configurations is also maximal 
for branches that propagate quasi-statically {v' 0), that is, when G' = G'g. Fig. [T| shows the 
dimensionless static energy release rate, 2nG'jKf = {Ff^ + F|i), as a function of the branching 
angle A. Note that this quantity equals 1/2 for "zero" branching angle and that it displays a 
maximum at a nonzero branching angle. 




Figure 7: The dimensionless static energy release rate immediately after branching 2fiG'g/Kf as 
a function of the branching angle A, when the primary crack is under mode I loading. 

The Griffith's criterion together with the principle of local symmetry (as given by Eq. ([8])) allow 
to determine direct dynamical properties of the branching instability. To be precise, Eq. (|18|) 
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reduces to 



m 

nvc) 



gi{v,) = Fl^{K). 



(19) 



This equation selects a critical velocity for branching, denoted by v^- When the fracture energy 
is velocity independent, Eq. (fT9|) simplifies to gi{vc) = Ffi(Ac) ^ 0.56. Obviously, the critical 
velocity depends on the material properties through n and cr only, but as shown in Fig. [HI it has 
a weak dependence on k, and can be reasonably taken to be Vc = 0.46c/j - its value for k ~ 3. 



0.46 



,0.45 



0.44 



0.43, 




Figure 8: The critical velocity for branching Vc as a function of k for constant fracture energy. 

When taking into account the velocity dependence of the fracture energy, Eq. ([8]) which se- 
lects the branching angle Ac is not modified, while Eq. (jl9p shows that th e critical veloc i ty for 



branching varies. In genera l, V{v) is an increasing function of the velocity (jBoudet et al 



1996 



Sharon and Fineberei . Il999l ) . Therefore, the left-hand side of Eq. (jl9p decreases faster than in 
the constant fracture energy case, and the energy balance can thus be achieved at a lower veloc- 
ity. Although T[v) can be a nonlinearly dependent function of the crack tip speed, it is only the 
amount of T{vc) /T{{)) which is of importance in determining Vc- In Fig. [9l the critical crack tip 
speed for branching is plotted for different values oiT(vr) and co mpared with experimental 
values for Glass and PMMA as given by (ISharon and Finebergl . 1 19991 ) . 



It is interesting to mention he re that if instead of the principle of local symmetry, the maximum 
energy release rate criterion (jErdogan and SihI . Il963l ) is used, then the equations determining 
the critical branching angle Ac and the critical branching velocity Vc are changed and so do the 
resulting values. In that case, Ac is given by the maximum of Ff^(A) + i^2i(^)' is just 

the solution of gi{vc) = F^^i(Ac) + F2i[Xc)- Interestingly, in practice Ac and Vc derived from 
these conditions correspond to almost the same values as those obtained from the principle of 
local symmetry. This is con sistent with other examples g iven in the literature, such as the case 
of kinked cracks treated in (jAmestov and Leblond ^ 19921). wh ere these two criteria yield almost 
the same numbers, and in contradiction to (|Adda-Bedial . l2005l ) where somewhat different results 
were obtained due to numerical inaccuracies in the determination of the functions Fim{X)- 



13 




r(v^)/r(0)-i 



Figure 9: The critical branching velocity for the case of velocity depende nt fracture energy and 



for AC = 3. The experimental values shown are estimates taken from (jSharon and Fineberg 



1999|) for Glass and PMMA. 



2.4 Dynamics of the branches 

The branching picture adopted here, i.e. with a fast moving main crack that stops, branches, 
and then the branches re-accelerate, might be questionable. Arguments criticizing this scenario 
rely on the fact that in experiments such full stops are not observed. However, this discrepancy 
can be related to different effects, such as the three dimensional nature of the experiments or 
to the fact that real materials are not ideally brittle, and so plasticity could smooth the present 
picture. Even though velocity fluctuations can be observed, if the crack stops completely then 
one would expect to see the branches accelerating gradually from zero velocity to Vc where the 
next branching event could take place. As we show below, the acceleration of the new branches 
can be rather large so that the rapid variation of the speed can be easily missed due to lack of 
temporal experimental resolution. 

In order to provide an estimate for the acceleration of the branches' tips, we use Griffith's energy 
criterion and the principle of local symmetry at each instant of the propagation of the tip of the 
branch. For a velocity independent fracture energy F, this yields 

1 



^ ^ , Gi (A,) (A.) - G2 jXc) Hu (Ae) T 



(20) 



H21 (Ae) Fii (Ac) Ki 

where terms up to order ^/I in the expansion of the SIFs ^ and the expression for the curvature 
a, Eq. ([9]), were used. In addition, knowing that immediately after branching, that is when 
i ^ 0, the crack speed of the branches is vanishingly small, i.e. v' 0, Eq. (|20p yields 
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By equating the last two expressions we get 



51 K) 



^1 



(22) 



where the numerical coefficient is computed by using the values of Fim, Gi and Him at A = Ac- 
Eq. ()20p is valid for a small extension of the branch, thus the velocity v' is small and one can 
safely develop 

9i{v') 



1 



~ C{k)- 



(23) 



51 (^^0 CR 

where C(k) is a numerical coef ficient o f order unity that can be computed from the asymptotic 
analysis of gi{v) as given in (jpreund . 1990 ). For a typical material value k = 3, one has 
C(3) = 1.15. Putting all together, one gets 



— ~ -0.6—crVI 
at Ki 



(24) 



which is a differential equation for £(t) that can ea sily be solved. The quantity T/K^ has been 
computed in Sec. for t he experimental setup of (jFineberg et al.l . Il992l : ISharon and Finebergj . 

2003 ). Using the result T/Ki = -0A3/VW for k = 3, the length of the 
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Livne et al 



branch extension is then given by 



The branch velocity is then just 



R+2 



i{t) ~ 0.017-^t 



(25) 



and finally the acceleration is 



v'{t) 



dt 



0.034-^t 
W 



dv' , Cp 

— ~ 0.034-^ 
dt W 



(26) 



(27) 



It is now obvious that the estimated acceleration of the crack tip after branching is very large 
{cr ~ 10^ — lO^m/s and W ~ 10^^ — lO^^m). This result might explain why even if the crack 
speed is vanishingly small immediately after the branching event, it would be difficult to detect 
it. 



3 Resolution of the static branching problem 



In this section we describe the analytical method we used for the resolution of the elastostatic 
problem of a long crack with two side branches. First the problem of straight branches is 
solved, rendering the universal functions and Gi. Then, the problem of curved branches 
is addressed, resulting in Hi^- The numerical results of this section were already summarized 
above in Figures [2l [3] and HI 
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3.1 Straight branches configuration 



Let us start by giving the general solution of the elastostatic problem depicted in Fig. [TOk . 
An infinite sheet is stretched in the presence of a crack contour consisting of a main crack of 
length L and two symmetric side branches of equal lengths £ emerging from a common origin. 
The angle between the two side-branches is denoted by 2\tt with < A < 1. In particular, the 
case of a main crack with two side-branches of infinitely small lengths is studied. The elastic 
potentials of the planar problem are determined for this geometry and loading. A conformal 
mapping of the exterior of this star shaped crack into the exterior of a unit circle allows to 
obtain integral equations for these potentials. Expressions for the stress intensity factors are 
derived. In the following, a detailed resolution of the mixed mode I-II loading is presented. 
Ac tually, the approach is analog ous to the kinked crack problem which was studied previously 



in (Amestov and Leblond 



19921 ^. 
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Z-plane 
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Figure 10: Conformal mapping of a star shaped crack in the Z-plane onto the exterior of the 
unit circle in the z-plane. 



3.1.1 Conformal Mapping, Potentials 



According to Muskhelishvili (jMuskhelishvilil . Il953l ) , the stresses and displacements at a point 
Z = X + lY = uj{z) (i.e. a two dimensional region - see Fig. [TO]l can be expressed, in the 
z-plane, in terms of the elastic potentials ^{z) and ^{z). We consider here the case where 
traction free boundary conditions are taken on the crack surfaces and where the loading is given 
by external stresses afi , and applied at infinity. The goal in this section is to solve for 
the Mushkelishvili potentials <I> and ^ outside the "star" shaped crack with two symmetrical 
straight branches, by using a conformal mapping transformation. The crack is located in the 
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Z = X + iY space (see Fig. [T0|) . it corresponds to the curve C, and its tips are located at 
the points = —L (tip of the left end of the original straight crack), B2 = £exp(iA7r) and 
Bi = iexp^—iXir) (tips of the branches). The region exterior to the crack is named . The 
potentials satisfy the following equation on the crack line C: 



$(Z) + Z<^'{Z) + ^{Z) = Const , (28) 

and the following boundary conditions at infinity: 

^{Z) = TZ, (29) 
l'(Z) = T'Z (30) 

with r = {a^ + (722)/4, r' = ((722 ~ <^ii)/2 + ^^12) given in terms of the stresses at infinity. 



The following conformal mapping (jSmithl . Il968l ) : 



Z = .{z) = ^{z- z,) {z - .3) f^)' (^)' (31) 

Z \Z-ZiJ \Z-Z3J 

maps the exterior of the "star" shaped crack described by the points Ai, Bi, A2, B2, A3, B3 in 
the Z-plane (see FIG. flUb ) to the exterior of a unit circle in the z-plane, with corresponding 
points zi,yi, Z2,y2, Z3,y3 (see FIG. [TOb). These special points are located at: zi = e~*", Z2 = 1, 
Z3 = e*", yi = e~*^, 2/2 = e*^ and ys = —1. The constants C, a and (3 are given in terms of the 
lengths L of the main crack, and i of each branch through the following equations: 

C = ^ [cos {a/2) f'-'^ (32) 

£ = -X{l-cosa)i X ^ (33) 

2 ^ ^lA(l + cosQ)i ^ ^ 

sin(/?/2) = yAsin(a/2) (34) 

In the z-plane, Eq. (f28ll for the elastic potentials reads: 

uj{z) 



(f){z) + —-^(j)'{z) + il){z) = Const , (35) 
u;'{z) 

{(j){z) = ^{Z), ipiz) = ^{Z)) and the boundary conditions at infinity become (j){z) = TCz, 
ipiz) = T'Cz. The quantity uj{z)/uj'{z) appearing in Eq. ([35]) takes different values in different 
sections of the crack line: 



{-1 + (1 - e-'^^nh{z) + (1 - e'^^^hiz)} Q{z) , (36) 



with Ii^2{z) = 1 if 2; belongs to Ci_2 respectively, and zero otherwise. Ci^2 are the branches of 
the crack, i.e. these curves unite the points Ai,A2 and A2,A3 respectively (see Fig. [TOb): and 
the function Q{z) is the following: 

Qiz) - (^-e-)(^-0(^-l) .37^ 
^" z(z + l)(z-e^/3)(z-e-/3)- 
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In order to solve equation (j35p . the following Lemma will be used: if / and g are complex 
functions defined and continuous in UC, analytic on O^, (including the point at infinity) and 
such that f{z) = g[z) for z E C, then / and g are constants and conjugate to each other. Thus, 
Eq. (j35p will be written in the previous form, i.e. as an equality between an analytic function 
and the complex conjugate of another analytic function. This is accomplished by defining: 



If z+ and represent points just inside and outside of the unit circle respectively, then by 
Plemelj's formula: 

XiM) - Xi,2{z-) = (1 - e^'^^-^)Q{z)W) (39) 

if z+j z~ are on Ci or C2 respectively, and zero otherwise. From equations (j37p . (jSSp one sees that 
convergence issues at the points e*'^ and e~*^ (poles of Q{z)) can be addressed by understanding 
these integrals with these poles slightly displaced into Vl~ (i.e. (5 ^ (3^iein the poles at e^*'^, 
respectively). Using these definitions and results in Eq. (j35p . the latter becomes: 



[z)-xi{z )-X2{z ) = Q^{z)(l)'{z)-xu{z-)-X2*{z-)-ij{z) + Const (40) 



where /*(z) = f{l/z) is an analytic function of z if f{z) is analytic. Thus, 

^*^>~ {z-e^P){z-e-^f^){z + l) ^ > 

Notice that if Q{z) has poles aX z = e'^e^^^ (i.e., in fl^), Q^,{z) has corresponding poles at 
z = e~''e^^^ (i-e., in Q'^). The condition of analyticity at 00 of the Lemma has to be examined 
for Eq. ([inD. Indeed the left hand side and the conjugate of the right hand side of the latter 
equation behave at infinity as: 

Hz) - xi{z) - X2{z) ~ TCz (42) 
Q.{z)(l)'{z)-xi*(.z)-X2*(.z)-i;{z) ~ -{T + T')Cz. (43) 

The behaviors at infinity are regularized (linear terms in z are eliminated there) if one adds the 
terms —TCz + (F + T')C / z and —TC/z + (F + T')Cz to the left hand side and right hand side 
of Eq. (j40p respectively, which then becomes: 

Hz) - xi{z-)-X2{z-)-TCz+{T + r)C/z 



Q.{z)(t>'{z) - xu{z-) - X2*{z-) - i'iz) - rC/z + (F + T')Cz + Const (44) 



Applying the Lemma to Eq. ([UD (left hand side), and using Eq. ([5^ . one gets: 

(F + r)c 



TCz 



z 



27rz 7ci {'n - z) 2711 Jc2 [rj - z) 
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i.e. an integral equation for (l)'{z) which can be written in the form: 

0'(z)=0'o(z)+£(0'(z)), (46) 

with 

^^(^)^rC+ ^^+P^ , (47) 

and the operator C is defined through: 



ru( (1 - e-^'^n f , (,y-e— )(7?-e^-)(??-l)/(r/) 



27ri 7ci ?7(ry + l)(r7 — e*^)(r/ — e~*'^)(r/ — z)2 

(l-e^ I ^-"X" -;-)(" + C»n»* (48) 

27ri ica ?y(T/ + l)(r?-e*/^)(?7-e-*^)(??-z)2 ^ ^ 



3.1.2 Expansion in powers of the crack extension length 



In this section we speciahze to the case i ^ 0, i.e. a situation with a long macroscopic crack 
with two micro-cracks right after branching. First, we write in the limit ^ — > an asymptotic 
expression for the constants C, a and /3 that follow from Eqs. ([32]) - (p^ : 



C 



a 







1 (l-A)L 





a 



(49) 

(50) 
(51) 



Writing 



and 



VLUiC)+aLV{C)+0{a'^ 



(52) 



equation (j46p becomes to order a (Notice that 17 corresponds to the lower half plane in the 
complex C plane, i.e. the poles at z = e^'^(l =F ^e) are now located at C = =fV^ — ie): 



+ aViO 



J2Xtt\ 



(l_e-i2A7r) ^^h{h'^ -l)(u{h)/VL + aV{h) 



ia—C + 
4 



1 



{h^-x){h-cr 



+ 



1 /i(/i2 - 1) (U(h)/VL + aV(h)) 
dh — ^ ^ 







which yields the following equations for U{C,) and V{C,)'- 

v{C) = v'^'HO+AViO 

where the functions C/*-^-' and V^^^ and the operator A are the following: 



L 



r + 



r\ {Ki-iK2) 



(53) 



(54) 
(55) 



(56) 
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|/I»>(C) = -.^c 



(57) 



"-^ 7-1 (/i2 - A)(/i - C)^ 47rz Jo {h^-X){h-Cy 



4m 



where 



Ki - iK2 = {a^^ - i<2) 




(59) 



are the stress intensity factors of the original single crack of length L under the same loading, 
and T is the non-singular stress (axx) at the original crack tip. 



Andersson's formula ( Andersson 



19691 ) for the stress intensity factors can be applied at the 



upper crack tip (one should get an analogous result at the lower tip), as follows: 



(60) 



where 6 is the angle between the X axis and the tangent to the crack at its tip, i.e. vrA in this 
case. To first order in q, we get 



A 



1-A 



AXn 



Also, writing the stress intensity factors up to order \/? as: 
one obtains: 



2v^e 



-iAti 



/i-aV/2 



— 



K 



(1/2) 



(1/2) 



1-A^^ 



(61) 

(62) 

(63) 
(64) 



Equations ([63]) and ^ show that the K^'s and A'/^^^^'s can be determined indepen- 

dently, i.e. the K^^^'^^^s depend only on the function U which can be found from Eqs. (I54p . (l56p 
and (j58p where the function V does not appear. Similarly, the k9'^'^^''s can be found from 



Eqs. ([55]) ■ ([57]) and ([58]) where the function U does not appear. This remarkable property holds 
only in the limit ^ ^ 0, and is easily evidenced thanks to the addition of the term e~*"^ in the 
expansion for <l)'{z) (Eq. ()52p ). Also, since the function C/^*^) depends on the three components of 
the stress tensor at infinity only through two parameters, namely the SIF's at the initial crack 
tip, the same holds for U and for the -ft^fs, i.e. they depend on Ki,K2 (this property is again true 
only in the limit £ ^ 0). In the same spirit, one could have thought that since depends on 
T and on a'^ the same is true for V , namely that it depends on bo th quantities T and a?^ . How- 
ever, this fact contradicts a universality prediction presented in (jAmestov and Leblondl . Il992l ) 
which states that the K^^^'^'^^s should only depend on the T-stress. In fact, the contradiction 
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is only apparent since it can be shown that the part of V which arises from (i.e. a V that 
solves Eqs. (j55p . (j57p with T = 0) has the following closed form 



[V (C)] 



T=0 



12 



c 



(Cen- 



(65) 



Now, it is easily seen that the function \y (C)]t=o T'Gxo at the point Q = \/A, so that o"]^ 
does not contribute to the SIFs K^^^'^'^^s which are given by expressio n (I64p . This observation 



shows that the universality predicted by ([Amestov and Leblondl . 11992 ) is respected. The exact 
result of Eq. (|65p serves as a useful check on the correctness of any numerical analysis of these 
equations. 



3.1.3 Numerical considerations 



Unfortunately, an analytical solution such as (jGSp was not found for the function U and for that 
part of V which is proportional to T. Therefore, in order to determine the functions Fim{\) 
and Gi{\) it is necessary to solve numerically for U and V . Let us begin with [/, where a useful 
decomposition is given by 

U{0 = -^(KiUiiC) - iK2U2{C)). (66) 
Vovr 



Eq. (j54p is now decomposed into two equations 

- - e"^^^^ r ri{rf -l)U^{ri)d'q 
Jc+ {if - A) (?7 - Cf 
.,2i\-K r riiif -l)Th^{r])d'q 



4i7r Jc+ iv^-X) {ri-Cf 



(67) 



where we deformed the integration paths away from the poles ±-v/A — ie onto two semi-circles 
denoted Cj^ and C2 respectively (i.e. |C ± ^| = ^, Im C > 0, oriented from —1 through 0, and 
from to 1 (see Fig. [TT]) ). 



Notice that the function f{z) = /(z) has been introduced, which is analytic when f{z) is analytic 
{f{z) coincides with f(z) on the real axis, and it is its analytic continuation into the rest of the 
complex plane). Eqs. (|67p can be solved numerically using an iterative method. Beginning with 
the non- homogenous term in the equations Ui^2 = 1, we can iterate using the operator A (given 
in eq. (j58p ) and obtain the following formal solutions 

00 

C/i(C) = 5:^"C/f^(C) (68) 

n=0 

00 

U2{c) = E(-ir-^"f^f (0 (69) 

n=0 
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Figure 11: The C-plane with some useful contours. 



which can be seen to converge, and to yield the required result. The convergence is due to the 
fact that the operator A is contracting in the space of functions defined and continuous on 
with, at most, weak singularities at 0, ±1 (see Appendix A for a proof). 



The simplest way is to compute the functions C^i(C) and U2{C) on the lower semi-circles 
and C2 respectively (see Fig. [TT]) since then C/i(C) and U2{C) which are needed in the integrals 
([67]) are trivially obtained on Ci U C2 by conjugation of Ui{() and U2{C)- In practice, we 
parametrized each of these two functions by two functions defined on [0, vr], namely 1^12(1) = 
Ui^2{^^^Y~^) on Cf and f/i_2(7) = ^1,2( ^^2'^'^ ) ^2 (with 7 e [0,7r]). In order to be more 
specific we write down explicitly the four singular integral equations for U12 that we solved 
numerically 



f/i%(7) = l ± 



1 - e 



-2iXiT 



8tt 



1+e- 



1+e- 



U12 {0)e dO 



± 



f^i,2(7) = l ± 



Svr 



U\2{0)e~'^de 



1 - e- 



1-e' 
2 



-A 

l+e- 



c 



Svr 7o (-i±f^j -A 



1-e'T 



(70) 



(71) 



Once ^7i(C) and U2{C) are known on the lower semi-circles, we use eqs. (j67p once again to obtain 
[/i(VA) and C/2(VA). Then, using definition (jj]) together with eqs. (|63p and (j66p we can extract 
the -Fjm's, namely the matrix elements relating the stress intensity factors immediately after 
branching at the tip B2 of the infinitely small side-branch to the stress intensity factors of the 
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main crack of length L in the absence of the side-branches at the leading order 

F„m = -l^il^X" s)\e-'>-u4-m 



-iX-K 



F22 (A) 



(72) 
(73) 
(74) 
(75) 



Results following from a numerical calculation of these fu nctions were presented above in Fig. 
[21 In that figure, we superimposed the results obtained in (jAmestov and Leblondl . Il992l ) for the 
case of a kinked crack (with the same angle) for the sake of comparison. 

The next stage is to get the G;'s. For that we need to solve eq. (|55p . Recall that we can dismiss 
the term since it will not contribute to the G;'s (due to the exact result we presented in eq. 
([65]) ). and so by denoting 



^(C) = ^v^(C) 



we then need to solve the following equation 



y(C) = iC + 



1 



-2iXiT 



rj (jf — 1) y (rj) dr] 
c+ (772 -A) [n-Qf 



+ 



rj [rf^ — l)V (rj) dr] 



(76) 



(77) 



This equation can be solved using the same iterative method as above, by taking the nonho- 
mogenous term in the equation V^^^ {(,) = iC, and iterating it using the operator A 



v[0 = Y,A^v^^^ (C). 

n=0 



(78) 



In order to solve this we apply exactly the same procedure as for f/i^2(C)) namely deforming 
the integration contours to Ci U C2 , and solving for V^(C) along the lower semi-circles Q G 
CiUC2- As before, we parameterize F(C) using two functions defined on [0,7r], namely (7) = 
y(^i^) on and ^''(7) = V{^^) on (with 7 G [0,7r]). This results in two singular 
integral equations for V""'^ that can be solved, and finally yield the required V{^/X). Then, using 
definition ^ and Eqs. (j64p and (j76p we can extract the following expressions for the G^'s 



2 V 1 - A V A 



1 - A 



(79) 
(80) 



Results emanating from a numerical resolution of V were presented a,bove i n Fig. [3l In that 
figure, we superimposed the results obtained in (lAmestoy and Leblondl . Il992l ) for the case of a 
kinked crack (with the same angle) for the sake of comparison. 
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3.2 Curved branched extensions 



We now consider the problem of the curved extensions defined in Fig. [12] with the aim to 
determine the Hi^^s (defined by Eq. dH)). The curved extension of one branch is described in 
terms of the coordinates li,2) where Yi is the axis parahel to the initial slope of the branch, that 
forms an angle vrA with the X axis {Y2 i s perpendicular to Y^). In this fram e of reference the 



shape of the extension can be written as (jLeblond 



1989 



Amestov and Leblondl . Il992l ) 



Yo = aY, 



3/2 



+ 



o (y,^) , 



(81) 



where a is a curvature parameter. A fictitious branch is defined as a straight line joining the 
beginning and the end of the branches, it forms an angle vrA with the X axis. The length of the 
fictitious branch is denoted i, and £ is the length of the curved extension. Expansions of i and 
A are: 



i 
A 



+ 



x + -^ri. 

TT 



(82) 
(83) 




Figure 12: Schematic representation of a straight crack with symmetrically branched curved 
extensions. The fictitious straight crack around which the perturbation expansion is performed 
is drawn on the lower branch. 

We will solve below for the Muskhelishvili potentials by a first order perturbative procedure, 
with the curvature a being a small parameter. This analysis will yield an expression for the 
SIFs of the curved branches as functions of the fictitious variables i and A. This expression will 
not yield directly the required -ff^m's defined in Eq. (jl]), since one needs to express the latter 
SIF's in terms of an expansion in \/l and as functions of A by simply changing variables from i, 
A to i, A. It is therefore necessary to r ewrite Eq. dH) in terms of I and A using Eqs. (I82p - ()83p . 



Following (jAmestoy and Leblondl . Il992l ) we get 



J2 Fim{X)Km+ [Gm{X)T5im + aHi^{X)Km]^l + 0{l) 



m=l,2 



m=l,2 



m=l,2 



m=l,2 



^^(A)^^^^ + aHUX)Km - -FUX)K„ 

vr 
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(where second order terms with respect to a have been disregarded). This expression is of the 
form 

(85) 



K'iii)= K[{1) +aJ2Him{X)K^\^i+0ii), 

a=0 ^ — ' 



where 



a=Q 



is the rth SIF at the tip of a straight extension of length i in the direction A, 



and the functions Him are defined by Him = Him 

HimW = Him{\) + 



F/^/vr, so we can invert it to get 



vr 



(86) 



The expression ()85p for the SIFs, which is exact to the first order with respect to a, is precisely 
of the form which will result from the perturbative analysis. It will therefore be easy to identify 
the functions Him, and the Him's will follow using Eq. (j86p . 



3.2.1 Perturbative analysis, integral equations 



The equations of the problem with curved extensions in the physical Z-plane take the same 
form as in the case of straight extensions, where again ^ and ^ denote the real Muskhelishvili 
potentials. We associate some new potentials and ^'"^ with <I> and ^' through analytic 
continuation by shifting the curved extension into the fictitious straight one. and '^^ have 
discontinuities across the fictitious straight extensions, while $ and ^ have discontinuities across 
the actual curved extensions. The values of and at each side of the fictitious straight 
extension are analytic continuations of the values taken by $ and ^ on both sides of the curved 
extensions. If Z = Zp represents the points of the straight fictitious extension (Z^ will represent 
points on opposite sides of the extensions), and if Zp+rju{Zp) represents the points of the curved 
extension of the upper (n) branch, then the original and shifted potentials are related through: 



$ (Z± + r,u{Zp)) ~ $^(Z±) + i^^{Zp)^"^{Z$) + O (rf) . 



(87) 



Also, and ^'"^ are expanded in powers of i]: 

= $o + ^i + 0(??2), 
= ^0 + ^-1 + 0(??2) . 

Combining Eqs. ([87D and (l88l) : 

$ (Z| + Vu{Zp)) ~ MZ$) + MZ$) + ^'o{Z$)Vu{Zp) 



(90) 



Following the perturbative analysis of (jAmestov and Leblondl . 11992 ) one obtains the following 
forms for Eq. (j28p written in z-space to orders 0(1) and 0{r]) respectively: 



Const 
Const' 



z + 



u;'{z) 



(91) 



uj'iz) 



<Poiz) , Mz) 
io'{z) u'{z) 



viz) 
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+ 



u;'{z) 



rj{z) 



(92) 



for z G C, with C the unit circle in z space corresponding to the fictitious star shaped crack 
{(po{z) = ^o{Z), etc.). Notice that Eq. is equivalent to Eq. ([35]), and that Eq. (01]) follows 
from the latter. Applying the previously mentioned Lemma (in section [3TT| right after Eq. (I37p ) 
to Eq. (I44p (to its right hand side, and considered for the fictitious crack), one obtains: 



where 



M^) = (r + T')c'z - rc'/z + Q4z)^'o{z) - xf^{z) - x£^(^) , 

Q*iv)4>oiv) 



(1 - e±2iA^) 



27ri ^ Jcia'^^ V{z-v) 



(93) 



(94) 



and with Q^{ri) given by Eq. (j41|) (with fictitious parameters). The poles of Q^{ri) at e*^ and 
e~*^ are displaced into 0^ (for Q{rj) they are at Vl^). 

The problem now is to solve Eg. (1921). Analys i s of t he potentials close to the singular points 
follows closely the one of (jAmestov and Leblondl . ll992l ). In order to apply the Lemma to Eq. ([92]) . 
one defines the functions: 



xji(^) 



Jci,2 - Z) 



2m 



i.e. by Plemelj's formula, they do satisfy: 

XiUz^)-X%{z-) = {l-e 



(1) (.+\ _ ..(1) (.-\ - n _ ^^■^i'^^)Q{z)(t)'^{z) , 



(95) 



(96) 



if z belongs to Ci^2 respectively, and zero otherwise. Also one defines: 

dt 



0?«(.) 



1 



2m Jci ^ {t - z) 



ri{t) 



u;{t)<p'iit) u;it)u;"{t)<p',{t) , Vo(i) 



+ 



vit). 



(97) 



Plemelj's formula can be applied in an analogous way to this equation. Applying these previous 
equations, Eq. ([92|) can be written as: 



ct,, {z)-x^^\z-)-x?{z-)-cp',^'\z-)-4^'\z-) 

= Q*{z)Mz) - xS{z~) - xtHz-) - 0?r^(z-) - 4*\z~) + Const. (98) 
for z G C. The Lemma implies that: 

Mz) - x'lhz-) - X^\z~) - 4>',^'\z-) - 4^'\z-) = Const. (99) 
Differentiating this previous equation, one obtains: 

cP[iz) = 4^'^'{z)+4'^'iz)+Wiiz), (100) 
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where C is the operator of Eq. (08]), with A ^ A. 

The steps in order to obtain the stress intensity factors at the tips of the curved extensions are 
the following: 

• i) Solve for </>o(2:) from the following equation: 

(t>',{z) = C\z) + C{<l^',{z)) , (101) 
(that follows from from Eq. (j9ip ) with 

{z) = rc7 + (r + r)c/z\ (102) 

(This step was practically done in section ISTT] with A and i instead of A and L) 

• ii) Evaluate ijjQ^z) from Eq. ([93]) . 

• iii) Calculate (i)%\z) from Eq. (f97]). 

• iv) Solve the integral equation (jlOOp for 4>'i{z). 

• v) Get the stress intensity factors using Andersson's formula (|60p with (j)' replaced by 
3.2.2 Expansion to order 1/2 in the extension length 

We will now expand the preceding equation up to order , or equivalently up to order a ( d is 
given by Eq. (|50p with A ^ A and £ ^ i). This will yield the functions Hi^- For that purpose 
we perform a change of variable, z = e*°^, and the following expansions of functions: 

0'o(z) = e-*°^[VZc/o(C) + a^^o(C) + 0(«')], (103) 
^'o{z) = e-*°^[VZXo(C) + 0(d)], (104) 
cj)[{z) = e-'^'^iVLUiiO + aLViiO + Oia^)]. (105) 

Expansion of the integral equation (jlOip for (jy^iz) up to order a leads to the following integral 
equations for Uq and Vq: 

UoiC) = U^'^HO+AUoiC), (106) 
^o(C) = vi''\c)+AVo(,(:). (107) 

where A is the operator A of Eq. §8^ with A ^ A, and C/^°^(C) = ^7d°\ ^o°^(C) = V^°^ are 
equivalent to those of Eqs. ([5^ - ([57]) (A — > A). Therefore, based on the results of section [3T] we 
can consider these as known functions from now on. 
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The expressions for Uq and Xq will also be needed here. The first one is obtained simply by 
differentiating Eq. (jl06p once with respect to C- In addition, we decompose Uo{() into Uoj{C) 
(j = 1,2) as is done in Eq. ([MD, i-e. f/o(C) = {KiUo,i{0 - i^2t/o,2(C))/\/8^, so we get for 

uLic) = ±^^^f^[ d,!i^l^MM^(l^f ..n^l^MM^ (108) 



2m Jc+ (r/2 - A)(7? - C)3 ^ni Jc+ (r?^ - A)(r/ - C)^ 

where the upper and lower cases on it correspond to j = 1,2 respectively, and convenient 
deformations of the paths r] £ [—1,0] and r] £ [0, 1] into ^}~^, respectively (in this way the poles 
at C = ±\/x — ie are avoided) as before - see Fig 111! 

In order to obtain Xq{C,) we differentiate and expand the expression for iPq{z) following from 
Eq. 



XoiC) = - + - + UoiC) - - u[^iC) 

V8^ 2(C2-A)2 "^^^ 2(C2-A) 

(1 _ e-2a.) ^0 ^^ hih^-i)u,(h) _ (l-e^-^-) h{h' - l)Uoih) ^^^^^ 



4vri J-i {h?-X)(h-CY ^vri Jo {h? - X){h - 

where now the poles are located at = iVT + ie (see the comment after Eq. (j4ip ). Therefore, 
in this expression the integrals over [—1,0] and [0,1] are deformed away from the poles into 
contours in the lower half-plane. However, using the lower semi-circles U as integration 
contours is not wise because then we encounter singularities in the integrand, as we are interested 
in evaluating Xq{C) on U C2 ■ It is wiser to use the two lower semi-ellipses defined by 
T] = (=f1 — (cos0 -|- isin^/2)) /2, < < vr (see Fig. [TT|) . In order to take advantage of the 
resolution of Uqj{C,) (j = 1, 2) done in the section [3TH we decompose Xq as in Eq. ([66]) 



Xo(C) = -^(i^iXo,i(C) - iK2Xo,2{0) . (110) 



And we get the following two expressions (j = 1, 2) 

. (C^ + f 1 - 3A)) + A CiC'^ - \) 

XoAC) = (-ir^ ^- Jr, UoAO-^f^K,iO 

2{e-X) 2(C^-A) 



-2iXn r r, {r,^ - I) Uo,j iv) dr, 1 - e^'^^ f V iv^ - ^) Uo,j (v) di] 



4i7r Jr- rf - \ {r/ - C)^ ^ivr Jr^ rf - \ [r] - 

We will now expand the integral equation (jlOOp for (j)'^ in powers of a. The first step is to get 
an expansion of io{z) and its derivatives: 

= -LaHe-Di-^YiArY • (112) 



^ f(C-.,(l)\l)\ 
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An expression for the gap rju{z) between the fictitious and curved upper branch is the following: 



Since 



it follows that 



with 



r]uiz) = -iae'^^\Z\{^Ji- J\Z 



\u;{z)\ = \Z\^^a^\e-l 



C 



(l-A)L/ A 



1-A. 



C + 1 

A/2 

a 



C-1 



C + 1 



C-1 



c 


A/2 


c 


A/2 


c + 1 




c-1 





1-A 



A/2 



1 - A. 



(115) 

(116) 
(117) 

(118) 
(119) 



All the necessary ingredients for expanding (/)2^^^ (z) and the integral equation (llOOl) are now 
available. Starting from Eq. (j97p . one gets the following expansion 

"1 hlh2 



0f)'(.) 



aL& /-i^^ h\h^-l\G{h) 



8tt Jo {h^-X){h-Cy 



{Uoih)+Xo{h) 



J2Xtt 



+ 



2(/i2 - A) 



(3/i^ - (1 + 5A)/i2 + 4A2 - A) 



ih' - A) 



Uo{h) 



(120) 



with a pole at ^ = v A — ie. Similarly, at the lower branch 



rii{z) = —iae 



r X 3/2 



e-1 



c + 1 



c-1 



(121) 



and then 

^o(i)' 



aha /-o h\h? - l\G{h) 



Svr J-i [K^-\){h-C)' 



[Uo{h)+Xo{h) 



i2Xi\ 



+ 



2{h? - A) 



h{h^ - l)f/^(/i) 



(3/1^ - (1 + 5A)/i2 + 4A2 - A) 



Uo{h) 



,(122) 



(/i2 - A) 

with a pole at = — In order to calculate more efficiently the integrals in Eqs. (I120p . 

()122p . we would need to deform the contours (as done many times before) into C^2- l^o^ that 
purpose the function G(C) may be continued analytically into disks that encircle the segments 
C G [— 1;0] and C G [0, 1]. These continuations are respectively: 



G{C) 
G{C) 



^_,,(l+A)/2^A(^ + l)(l-^)/2(C - l)(l-^)/2 - Vl - A 
g-i.(l-A)/2^A(^ ^ l)(l-^)/2(C _ l)(l-A)/2 _ y^ry 



A 



1 - A^ 
1 - A 



A/2 



A/2 



(123) 
(124) 
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The following step is to solve Eq. (jlOOp for (t)'i{z). Replacing the expansion for 4>'i{z) in Eq. (jlOSp 
into Eq. (jlOOp . one gets that (f)i2 (-2) are of order a. As a result the equation for Ui (analogous 
to Eqs. (|lU6p - (|lU7p for [/q and Vq) becomes Ui{C,) = ^f^i(C) (where here too A is the same 
operator as in Eq. (|58]) . but with A ^ A). Due to the contracting nature of A (see Appendix 
A) this implies that C/i = 0. Expanding Eq. (jlOOp to order d, using the previous expressions for 
(i)i^^'{z) (Eqs. ([I22|), (11201)), we get 

^i(C)=^/°^(C)+^^i(C), 



(125) 



with 



v<»'(C) 



- l|G(/i) 

r""(/i2-A)(/i-c)2 



+ 



2(/i2 - A) 



/i(/i^ - l)U'Q{h) 



[U^{h)+X^{h) 

(3/1^ - (1 + 5A)/i2 + 4A2 - A) 



{h^ - A) 



[/o(/i) 



h\h'^ - l\G{h) 
"■"(/i2 _ A)(/i-C)2 



{?7o(/l)+Xo(/l) 



+ 



2(/i2 _ A) 
By further decomposing Vi{C) into 

ViiC) 



h{h^ - l)U^{h) 



(Sh^ - (1 + 5A)/i2 + 4A2 - A) 



(87r)3/2 



(/.2 - A) 
{KiVi,i - iK2Vi,2) , 



(126) 



(127) 



we get: 



^ h\h^-l\G[h) 

1 (/l2-A)(/l-C)2 



{c/oj(/i)±Xo,j(/i) 



-i2A7r 



2(/i2 _ A) 



h{h'-l)U'{h) 



(3/1^ - (1 + 5A)/i2 + 4A2 - A) 



h\h^-l\G{h) 

(/l2 - A)(/l-C)2 



±- 



'2{h?-X) 



h{h'-l)U',^^{h) 



(/.2 - A) 
(3/1^ - (1 + 5A)/i2 + 4A2 



A) 



f^o,i(/i) 



(128) 



{h^ - A) 

(the upper and lower cases of ± correspond to j = 1, 2 respectively) where previous decomposi- 
tions for C/o(C) and Xq{C,) (given by Eqs. (l66l) . (|110p ) have been used. Then, Eq. (|125p leads to 
the following equations to be solved: 

y/?(C) + ^Vi,i(C), 

(129) 



^1,2(0 = <2^(C)-^Vi,2(C), 

which are solved as in the previous section by iterations as follows: 

00 

n=0 

00 

n=0 



(130) 
(131) 
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In order to apply Andersson's formula (j60p at the upper branch tip, one uses 



(5 ~ vrA H — a\ 
2 



^(ei/3)(ei/3 + i)(ei/3_g-i/3) 



Also, since f3 ~ \/J,a, and A = A + aVi/ir we get 
Then, for the curved crack case one can write (Eq. (jSOp ): 



1-A 



[/ VA + VldV VA 



1 - A 



A/2 



^-inX^-ia\/~i/4 



A/2 



ttA 
a=0 



+2V^e^ 



- iAvr 



i-A 



1-A 



1-A 



A/2 



V, VA 



(132) 



(133) 



(134) 



(135) 



Using the decompositions given by Eqs. ([66]) . (jllOp and (jl27p . the last expression for the stress 
intensity factors becomes: 



K[ii)-tK!,ii)= K[{i)-iK'^{i) 



a=0 



-if/o,,|VA)+-.' ' ''"^ 



A/2 



vr V 1 - A \ A 



Via VA 



i^2 



A/2 



(136) 



The HimS are extracted from the defining expression Eq. ()85p together with the previous result. 
Notice that extracting the Him's from the last equation i nvolves evaluating the functio ns Vij{C) 
at \/x (on the real axis). This involves, as explained in ( Amestov and Leblond . 1992 ). crossing 
of the pole at r/ = \/x for V^/j^(C)- Appropriate account of this difficulty results in that if one 

■"^^ ' , , . fj.Qj^ ^j^g f/'g should be doubled. Thus, 



just replaces ( 
we get: 



in expression (|128p the contribution 

A/2 



5R<^ e 



zAtt 



A/2 



1 /i-a' 



A/2 



Vi,2 V A 



^21 (A) 

^22 (A) 



1 /l-A 



A/2 



4^/2 V A 



9< e 



-zAtt 



2f?7o,i 



A/2 



1 /^l-M cjl -.A. 



4v^ V A 



1 /i-a' 



A/2 



V^l,2 V A 



40) 
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Finally, Eq. (|86]) is used in order to extract the Him^s from the HimS. Results following from 
a numerical resolution of the problem were presented above in Fi g. HI In that figure, we 
superimposed the results obtained in (jAmestov and Leblondl . Il992l ) fo r the case of a kinked 
crack (with the same angle) for the sake of comparison (actually in (jAmestoy and Leblond . 
1992) only results up to A = 80" are presented, so we extended the results obtained there to the 



whole range A G [0, 1] in order to allow for a comparison between the two cases) 



4 Elastostatic analysis of a crack in a strip geometry 

(u =0 or a =0) and u =5 

^ X xy ' y 

^ y=1 

, u =0 

) „ y y=0 

a =0 ' 

xy 



Figure 13: A semi-infinite crack in a strip of unit half-width. The top and bottom edges of the 
strip are fixed at height ±5 above their initial positions, and either these edges are clamped 
(first problem) or free to slide (second problem). 

In this section we calculate the ratio T /Ki for a static semi-infinite crack in a strip of unit half- 
width in an isotropic elastic medium described by Lame coefficientes A, /i (see Fig. [T3|) . The 
system is loaded by imposing normal displacements ±5 to the top and bottom edges of the strip. 
We will consider both clamped {ux = 0, first problem) and shear free {axy = 0, second problem) 
boundary conditions. In both cases, the symmetry of the problem allows to focus on the upper 
half plane only. Also, we subtract the solution of the unbroken strip, which corresponds to a 
state of uniform stress. This contribution will be added when needed. The boundary conditions 
corresponding to the first problem are 

Ux{x,l) = Uy{x,l) = axy{x,Q) =Q , (141) 
ayy{x,^) = -cjoQ, < , (142) 

Uy{x,Q) =0, |x| > 0, (143) 

while the second problem has the following boundary conditions 

axy{x,l) = Uy{x,l) = axy{x,Qi) ={) , (144) 
ayy{x,f)) = -cJoo , kl < , (145) 
Uy{x,Q) = 0, |rE|>0. (146) 

Here (Too can be determined from the solution of the unbroken strip and its value is different for 
each problem. Notice that for the second problem, any addition of a constant stress (Txx{x — > 
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ibcxD, y) = Too is consistent with the boundary conditions, and thus with the formulation of this 
problem. In both cases we have 

Uy{x ^ -oo,y ^ 0'^) = 5 , (147) 

Introducing Fourier transforms in the x-coordinates as y{x) = ^ Y{k)e^'^^^dk, the elasto- 
static equilibrium equations can be solved without difficulties. Moreover, applying the bound- 
ary conditions p4ip and (|144p of the first and second problem respectively, allows to reduce 
the determination of the stress field to the problem of solving the following equation with the 
appropriate boundary conditions: 

(A + 2^)F(A:)^-^, (148) 

where 

rca 

T.yy{k) = ^yy{k,0)= ayy {x , O) c'"'' dk , (l49) 



tkx , 

Uyy{X,V)(i 

— OO 

oo 

Uyik) = Uy{k,0)=l Uyy{x,0)e'^''dk, (150) 



and F{k) is a known function that is obtained in a closed form for each problem. As will be 
explicitly shown below, the behavior of F{k) at large and small A;'s is given by 

F{k) ~ /o, |A:|«1, (151) 
F{k) ~ foo\k\, |A;|»1. (152) 

The piecewise boundary conditions (I142p . (ll43p and ()145p . ()146p suggest the use of the Wiener- 
Hopf decomposition method. Let us introduce 



oo 



^yy(^) = ayy{x,0)e"'''dx , (153) 

^yyik) = f ayy{x,0)e'^''dx , (154) 

J —oo 

with Uy{k) and Uy{k) defined similarly. Notice that T,~y(k) (Uy{k)) is analytic for < 
and T,+y{k) {U+{k)) is analytic for Qfc > 0. Using ([112]) . (fTi3l) and (fTi5l) . (fTi6l) together with 
Eq. (jl48p one thus obtains 

- (A + 2i,)F{k)Uy {k) = ^+y{k) - ^ . (155) 
Let us suppose that one can write 

= (156) 

where F~[k) has neither zeros nor poles for 9 A; < and F^{k) has neither zeros nor poles for 
Qk > 0. Then Eq. (jl55p is rewritten as 

-{X + 2fi)kF-{k)Uyik)=ki:+yik)F+ik)+iaooF+{k). (157) 
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The left-hand side of Eq. (jl57p is analytic in the lower half plane, while its right-hand side is 
analytic in the upper half plane, and both sides coincide on the real axis. By the theorem of 
analytic continuation, and in order to retrieve the square root behavior of the stress field at 
the crack tip, both sides must equal a constant. This constant can be fixed by examining the 
behavior of the expression for k ^ 0. In fact Eq. (jl47p implies that 

6 



U-{k)^ 



ik 



k^O . 



Equations (11570 and (11580 then imply 



Tlyy{k) 



yy 



ik 



' ik F+{k) ' 



where 



o-o = (A + 2fx)5 . 

Now, since we are interested in the T-stress we look at the difference 



1 

0+) - ayy{x,0+) = — G{k)T.yy{k)e 



ikx 



dk . 



(158) 



(159) 



(160) 



(161) 



where G{h) = — 1 + Ti^xik) /T,yy{k) will be simply a given function for each problem. We will 
also need the stress intensity factor (SIF) Ki in each case, for which the expression is given by 



lim 



V-2zkayy{k) = croV2/o/( 



(162) 



Before examining each problem separately, the method of decomposition F(k) as F~ {k) / F~^ (k) 
which is necessary for the Wiener-Hopf method will be presented. To do this, let us define the 
function 

H{k) = , ^^^^ . (163) 

where /o and foo are chosen to get the two limits \k\ 0, itcxD right for real k. The function H{k) 
is a bounded, even function, which ten ds to 1 for |fc | 0, ±00. Therefore, we can approximate 
this function u sing the method of Fade jBakeii . Il975l ). This method was first used in the context 



of fracture by (jBouchbinder et al 



2002 



A Fade approximant, is that rational function of a specified order whose power series expansion 
agrees with a given power series to the highest possible order. In the present case, the Fade 
approximation of g{k) will be of the form 



H{k) 



P2N{k) 
Q2N{k) ' 



(164) 



where P2N{k) and Q2N{k) are polynomials, which must be even, and 2N is the order of the 
approximation (in practice we use 2N = 30 which gives accuracy better than 10~^ for the 
desired quantities). Then, we find all the complex roots £1, ...,^27V and Ai, A2Ar of P2N{k) and 
Q2N{k) respectively (actually it is even a simpler problem to find the N roots of the polynomials 
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P2NiV^) and Q2n{V^), denoted by ri, ...,rjv and pi, ...,pN respectively, from which we recover 
i2n-i,^2n = and A2n-ii ^2n = ^^/Prl)■, and SO we get the factorization 



H{k) ~ 



1_ 



fc2 



(165) 



Pjv • 



The Wiener-Hopf decomposition of a function of this form may be carried out by inspection: 



F+{k) 




(166) 



This approximation schemes converges when taking larger and larger A^'s, and so the choice 
of is a matter of the desired accuracy. The a dvantage of this appro ach over expanding in 
Chebyshev polynomials, as done for example in (jLiu and Marderl . Il99ll ). is twofold. First, by 
using the same number of series coefficients a better approximation for g{k) is obtained. Second, 
when factorizing the expansion in Chebyshev polynomials to order 2N one has to find roots of a 
polynomial of order 2N while for the comparable Fade approximation one has to factorize two 
polynomials of order which is simpler. 



4.1 Solution of the first problem 

For the problem defined by the boundary conditions (|141|) - (jl43p . the function F{k) is defined 
by 

2 k [(^2 + 1) + 2(k - Ifk^ + [k' - 1) cosh(2A;)] 
^ ' k2 (k + 1) sinh(2A;) - 2A;(k - 1) ' ^ ' 

so that /o = 1 and /oo = 2(k — 1)/k^. Also, the function G{k) is given by 
Gik) = -4 (k- I)2fc2 + K 

^ ' (^2 + 1) + 2(k - 1)2A:2 + (k2 - I) cosh(2A:) ' ^ ^ 

By solving the problem of an unbroken strip with the same boundary conditions one gets 

f^oo = (To = (A + 2^J)5 . (169) 

Now, using the decomposition given by Eqs. (jl59p and (jl66p . one has 

N /l 

The stress intensity factor given by Eq. (|162p yields 



Ki = 2o-o\/k- 1/k ■ (171) 

The T-stress is related to the asymptotic value of {axx — (^yy) at y = for a; ^ with adding 
the solution of the unbroken strip 

2u 1 2 

T = lim [a,,(x,0+) - (7yy{x,0+)] - —^^0 = / G{k)^yy{k)dk - -ao . (172) 

In Fig. [T^ the ratio T/Ki is given as a function of k. Notice that T is always negative for this 
problem. 
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Figure 14: The dimensionless ratio T/Ki as a function of n for the clamped boundary condition. 
Notice that for k = 3, T/Ki = —0.43, or in dimensioned quantities T/Ki = — O.S/VTy, where 
W is the width of the sample. 



4.2 Solution of the second problem 



For the problem defined by the boundary conditions (|144|) - (jl46p . the function F{k) is 

, . _ (K-l)fc(2fc + sinh(2fc)) 

so that /o = 4(k; — 1)/k^ and /oo = 2(k — 1)/ . Moreover, one has 



(173) 



G{k) 



-4k 



and 



2k + sinh(2A;) ' 



1 + — ^ — 

K V do 



(174) 



(175) 



where Too = (Txx(oo). Now, using again the decomposition given by Equations (|159p and (|166p . 
one gets 

y/ i — IK/ A 

n=l 



Yjyyik) 



ik 



N I 



ik 



(176) 



In this case, the stress intensity factor is given by 

4(k - 1) 



-do 



and 



27r 



G{k)Y,yy{k)dk o"o 



2(k - 1) Too 
«; do 



(177) 



(178) 



Fig. [T5k shows the behavior of T/Ki as a function of k for various values of Too- For every value 
of K there exists a critical value of Too, denoted by Tcr, at which the value of T changes sign 
(see Fig. [T5b). 
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Since a change in the sign of T imphes a transition from a stable crack growth (T < 0) to an 
unstable one (T > 0), a destabilization/stabilization of the growth process may be induced by 
tuning the value of Too- Interestingly, this feature is not shared by the first problem (clamped 
edges) and is a particular property of the shear free boundary conditions. 
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Figure 15: (a) The dimensionless ratio T/Ki as a function of k for the shear free boundary 
condition. The different curves correspond to different values of T^/ctq. (b) T(.^/(Jq for various 
values of k. 
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Appendix A: Proof that A is contracting 



Take the space of functions which are defined and continuous on such that the norm defined 
by 



Max 



(c' - i) Cf (C) 



(179) 



is finite. Here, we intend to prove that the operator A (given by Eq. (|58p ) is contracting in this 
space, i.e. that there exists a constant c smaller than 1 such that 



P/ll<c||/|| 

for every function in that space. Rewriting A by shifting the contours to '^^ S^t 



(180) 
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Now using 
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5— 3 COS g S 
2 2 



and 



5+3 COS e g- I ^ gg|. following estimate 



P/ll < 11/11 -Max 



sinvrA 



CeC-„ L 47r 



(c^ - i) c 



d0 



[(-^)-c] (-^) - (^) 



+ 



.(182) 



We focus on the left branch Ci , where C 



l+e'T 



(7 G [0,7r]). There we get |(C^ - 1)C| 



5+3 cos 7 gjj^^_ £qj, £]^g integral in the last inequality, we are not able to bound it analytically, 
but it is not difficult to show numerically that it obtains its maximum at A ^ 1 and 7 ^ vr 
where it behaves like -^ji^:^- And so for every A < 1 we get that Max{- • •} < 1. Similarly, for 



CGC- 



the right branch C^, where C = 2 ^ ■^D exactly the same bound, since the only 

difference with respect to is 7 — > vr — 7. Therefore we get what we need, i.e. there exists a 
constant c < 1 such that Eq. ()180p is obeyed. 
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